Shear thickening in molecular liquids characterized by inverse melting 
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We studied the rheological behavior of a molecular solution composed of a-cyclodextrin, water 
and 4-methylpyridine, a liquid known to undergo inverse melting, at different temperatures and con- 
centrations. The system shows a marked non-Newtonian behavior, exhibiting the typical signature 
of shear thickening. Specifically, a transition is observed from a Newtonian to a shear thickening 
regime at a critical shear rate 7c. The value of this critical shear rate as a function of T follows an 
Arrhenius behavior 7c(r) = B exp{Ea/KBT), with an activation energy Ea close to the value of 
the hydrogen bond energy of the O-H group of the a-CD molecules. We argue that the increase of 
viscosity vs shear rate (shear thickening transition) is due to the formation of hydrogen bonded ag- 
gregates induced by the applied shear field. Finally, we speculate on the possible interplay between 
the non-Newtonian rheology and the inverse melting behavior, proposing a single mechanism to be 
at the origin of both phenomena. 

PACS numbers: 61.25.Em, 62.10.+S, 66.20.Ej 



The study of complex fluids, like colloidal suspen- 
sions, surfactant solutions, polymeric systems, proteins, 
glasses and granular materials has gained during the last 
decades, the increasing attention of scientists [liQ- Due 
to the varieties of interactions that give rise to a rich phe- 
nomenology, these systems represent the prototypes to 
understand most of the relevant problems in the physics 
of soft matter Q . Besides the wide interest for fundamen- 
tal research, the study of complex fluids is of great rele- 
vance in for practical applications due to the possibility 
of tailoring different system properties as for example the 
mechanical or thermal properties, or the capacity to solu- 
bilize and transport materials, etc. This class of fluids at 
variance with simple liquids is characterized, at molecular 
level, by the presence of complex micro-structures which 
give rise to non-Newtonian behaviors. Consequently the 
study of the rheological response of these systems, pro- 
vides critical information for their performances in many 
industrial applications. In fact, due to several details of 
the system like inter-particle interactions, volume frac- 
tion, electrostatic forces, steric repulsion, size and shape 
of particles, Brownian motion, hydrodynamic forces etc., 
complex fluids may or may not show the appearance of 
a dynamical instability. In this context rheological tech- 
niques represent a widespread method to monitor, inves- 
tigate and/or predict anomalous behaviors of complex 
fluids. Common rheological measurements allow to in- 
vestigate phenomena like shear thickening, shear thin- 
ning and shear banding. In general for a Newtonian fluid 
the viscosity is shear rate independent while for a non- 
Newtonian fluid it depends on the applied shear. For the 
latter it is intuitive to think that an applied shear could 
destroy existent fluid structures giving rise to a decrease 
of viscosity (shear thinning), however there are some flu- 
ids which behave in a completely different and counterin- 



tuitive way since when subjected to shear they manifest 
an increase of viscosity giving rise to the phenomenon of 
shear thickening. This behavior has stimulated the inter- 
est around this phenomenology. Many complex fluids, in- 
cluding dense colloidal suspensions and dilute surfactant 
aqueous solutions, exhibit a shear thickening transition. 
The continuous development of new theories and the in- 
troduction of modern concepts regarding shear thicken- 
ing and the physics of complex fluids and soft matter 
promote the realization of targeted experiments and the 
implementation of experimental techniques aimed to give 
complementary results. Recently, on the basis of the 
new knowledge acquired in the field of soft matter and 
glass transition a class of models for shear thickening and 
"jamming" in colloidal suspensions has been proposed [3| 
where "jamming" induced by shear indicates an arrested 
phase. In these models starting from the mode coupling 
theory (MCT) of the glass transition, flow curves and cor- 
relation functions arc calculated using memory functions 
which depend on density, shear stress and shear rate. In 
this connection in very recent viscosity studies 0, 0, 0| 
on suspensions of non-Brownian particles subjected to 
shear there is evidence of shear thickening and/or per- 
manent "jamming" . In a further recent spin model for 
shear thickening reported by Sellitto and Kurchan [8] a 
generic mechanism for shear thickening analogous to en- 
tropy driven phase reentrance is proposed. In the latter 
paper the analogy between entropy-driven transitions in 
which systems freeze upon heating and the transitions 
in which they jam under the action of stirring, is ex- 
ploited. Starting from an Hamiltonian modelling a liq- 
uid that undergoes inverse melting, they introduce a force 
field representing stirring and obtain a direct derivation 
of shear thickening and glass transition found through 
MCT by Holmes et al. for an inverse freezing sys- 
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FIG. 1: Phase diagram of aCD-water-4MP solutions from 
ref TgI] . The transition temperatures are plotted as a function 
of the concentration of qCD in 4MP. The arrows indicate the 
concentrations and the high temperature intervals at which 
the rheological measurements have been performed. 



tem. This phenomenon happens when a liquid heated 
at constant pressure undergoes a reversible liquid-solid 
transition, thus generating a solid with entropy higher 
than its liquid counterpart [l^, 11 , . 

The theoretical possibility to correlate two apparently 
independent system properties, as "shear thickening" and 
"inverse melting" , is exciting and promising, however a 
clear experimental evidences of the coexistence of these 
two behaviors in a single system was -up to now- absent. 
In this paper we present the first experimental evidence 
of shear thickening in a molecular liq uid that is known 
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to undergo inverse melting [13|, [IJ, 

We performed rheological measurements as a func- 
tion of temperature and concentration on molecular so- 
lutions of a-cyclodextrin (aCD) (Csg^feoOao), water and 
4-methyl-piridyne (4MP) {CqHjN). The solutions under 
applied shear exhibit the phenomenon of shear thicken- 
ing. In particular we show that the system goes trough 
a shear thickening transition at a critical shear rate 7c 
which increases with increasing concentration. It defines 
two different regimes: at low shear rates (7 < 7c ) the 
fluid is Newtonian-like, while at larger shear rates, 7 > 
7c, a transition from the Newtonian behavior to the shear 
thickening regime is observed. Viscosity measurements 
have been performed at the different concentrations and 
temperatures indicated by the arrows in the phase di- 
agram of fig. [T] with the aim to understand a possible 
connection between inverse melting and shear thickening 
and the nature of the transition between the two different 



rheological behavior observed. 

The sample was prepared by dispersing aCD hydrate 
in liquid 4MP 98% purity (both purchased from Aldrich) 
and dcionized water and stirring for about 4 hours un- 
til the suspensions had cleared. The measurements were 
performed with a Reostress RS150 rheometer (Haake) 
in Couette geometry equipped with two coaxial stainless 
steel cylinders. The inner one, the sensor, fixed to a ro- 
tor motor has a diameter of 2 mm (model Z20DIN). The 
sample was placed in the 0.85-mm gap between the two 
cylinders; the temperature of the outer one was stabi- 
lized through a liquid fiux cryostat (DC50-K75; Haake). 
The rheometer was used in control rate (CR) mode to 
measure the shear stress cr as a function of shear rate 7, 
and hence to determine the shear viscosity ?7(7) through 
the relation a(7) = ?7(7)7. The measurements were per- 
formed as a function of the shear rate 7 in the range 
0-^1400 s^^. The sample was prepared at nine different 
concentrations of aCD, water and 4MP with molar ra- 
tios l:6:x, respectively, with 50 < a; < 200 as shown in 
the diagram phase T — c of fig. [1] 16] . The measurement 
were performed in the temperature range 259-365 K. 

As examples, in fig. [2] the (apparent) shear viscosity 
77(7) = o'(7)/7 is shown as a function of the shear rate 
for few selected temperatures and three sample concen- 
trations: c = 125 mg/ml, c = 201 mg/ml and c = 143 
mg/ml. As exemplified in fig. [H the measurements are 
characterized by two different regions separated by a crit- 
ical shear rate value 7c: at low shear rates, below 7c, 
the solutions display Newtonian behavior (the slight in- 
crease is related to spurious inertia effects of the rheome- 
ter as ascertained by performing similar measurements 
on a "blank" -water- sample) while at high shear rate, 
above 7c, a transition to a clear shear thickening regime 
is observed. 

From each critical shear rate, 7c, a critical shear stress 
value, CTc, has been obtained as ctc = r?(7c)7c- Interest- 
ingly, as shown in the inset of fig. [3l a log-log plot of 
all the temperatures and concentrations values of CTc vs 
the corresponding 7c displays a single power law behav- 
ior (Tc ^ 7c (the value of the exponent turns out to be 
a = 2.28 ± 0.2). For each T and c point, the critical val- 
ues of viscosity (or the stress) and shear rate were used 
to normalize the viscosity curves reported in fig. [51 The 
normalized shear stress, (j/<Jc, as a function of normal- 
ized shear rate, 7/7C, are shown in fig. [3] for selected T 
and c points. In fig. [Hta) the curves at constant T (T = 
293 K) at different concentrations are shown. The data 
collapse on a unique master curve and, for 7 > 7c, they 



show a common power law behaviour |l8l . [l3 . 20l . 
in which the shear stress is expressed as: 
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(7(7) 



(1) 



where n ^ 1.8. The fit performed according to the model 
of Eq. [1] is plotted (full line) for 7 > 7c in fig. [3] and it is 
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FIG. 2: Shear viscosity versus shear rate of aCD-water-4MP 
solutions for three different concentrations of qCD in 4MP (c 
= 125 mg/ml, c = 143 mg/ml and c = 201 mg/ml) at the 
indicated temperatures. 



compared with the Newtonian flow corresponding to n=l 
and reported as dashed hne for reference. The same be- 
havior of the normahzed shear stress vs normaUzed shear 
rate is shown at fix concentration (c = 125 mg/ml) and 
different temperatures in fig. El^b). This common trend 
imply that all the measurements scale on a unique master 
curve for all the temperatures and concentrations inves- 
tigated determining the existence of a single rheological 
behavior for the fluid. 

In fig. m the jc — T dynamical phase diagram is shown, 
it is evident that the onset of the shear thickening occurs 
at larger shear rates as concentration is increased: the 
more concentrated the solution, the higher the shear has 
to be to observe shear thickening. Moreover the criti- 
cal shear rate decreases as the temperature is increased 
suggesting that the shear thickening transition is favored 
at higher temperatures. The data follow an Arrhenius 
behaviour: 



7c (T) ^ Be^ 



(2) 



where Ea and Kb are respectively the activation energy 
and the Boltzmann constant. This behavior seams to re- 
sembles the one found in solutions of micelles Bin, Hi, 



FIG. 3: a) Log-log plot of the stress versus shear rate normal- 
ized to their critical values (defined in the text) for different 
concentrations at T = 293 K. The power law fits of experimen- 
tal data are shown before and after the critical point (solid 
line). The Newtonian behaviour (dashed line) is also plotted 
for comparison, b) Log-log plot of the stress versus shear rate 
normalized to their critical values for different temperatures 
at c = 125 mg/ml. In the inset a log-log plot of the criti- 
cal stress as a function of critical shear rate is shown at the 
indicated temperature for different concentrations of aCD in 
4MP together with a power law fit (dashed line). 



but at variance with them a decrease of shear rate with 
temperature is observed. This apparently counterintu- 
itive phenomenology can be rationalized by noticing that 
in an inverse melting system there is a slowing down of 
the dynamics at increasing temperature. Let us con- 
sider a characteristic time dictated by the shear fiow, 
i. e. Tj = l/'jc and the characteristic relaxation time of 
the fluid, Ta, we should expect a "transition" when these 
two timescales start to "interact" . Thus, a decrease of % 
implies an increase of t^, consequently of Tq, and hence 
a slower dynamics induced by temperature. The result- 
ing fits with an Arrhenius law (dash-dot lines in fig. |4]) 
are in good agreement with the experimental data. The 
parameter Ea obtained from data interpolation is plot- 
ted as a function of concentration in fig. [51 It increases 
linearly with concentration suggesting that a greater acti- 
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FIG. 4: Critical shear rate of Q:CD-water-4MP solutions as a 
function of temperatures for some of the investigated concen- 
trations (symbols) . The Arrhenius fit of the data (as discussed 
in the text) is shown for some of the investigated concentra- 
tions (dash-dot lines). The regions where the solutions exhibit 
two different rehological behaviours are displayed. 



vation energy is required to activate the shear thickening 
regime at higher concentrations. In particular we would 
like to emphasize that the values found for the activa- 
tion energy, ranging between 0.1 and 0.2 eV, are of the 
same order of magnitude of the intramolecular hydro- 
gen bonds (HBs) energy, Ehb, of the aCD molecular O- 
H group (£^^5=4.5-4.9 kcal/mol = 0.196-0.213 eV) ^ 
and of water {Ehb ~ 5 kcal/mol) [27|]. This result seems 
to corroborate the hypothesis proposed in ref [l^ which 
suggests that the applied shear rate favours the brak- 
ing of intramolecular HBs at a critical value 7c, which 
depends either on concentration and temperature. This 
in turn could facilitate the reformation of intermolecular 
HB with water and aCD giving rise to the formation of 
more complex aggregations and to an increase of viscos- 
ity. This hypothesis is also reinforced by the temperature 
behavior of the critical shear rate of fig. 2] which decreases 
at increasing temperature. In fact, the increasing of T 
favours the thermal motion of molecules and the breaking 
of HBs implying that a smaller 7c is required to activate 
the shear thickening regime. In addition fig. |4] suggests 
that at fixed temperature for higher concentrations of 
aCD a higher energy is required to break the intramolec- 
ular HBs. Furthermore fig. [5] underline the role played 
by 4MP on the shear thickening transition and on the 
breaking of HBs, it shows in fact, that the activation en- 
ergy required for the transition decreases at decreasing 
concentration of aCD and departs from the value of the 
Ehb for pure aCD {Ehb ~ 0.2eT^). 



FIG. 5: Activation energy obtained from the fit reported in 
fig|4]as a function of the concentration (symbols). A linear fit 
of the data is also shown (dashed line) . 



These results can also be interpreted considering a time 
scale r determined by the inverse of the critical shear rate 

= l/7c. When 7cT < 1 the aCD molecules subjected 
to shear have enough time after a collision to move with- 
out remain H-bonded to other molecules. At increasing 
shear rate the condition 7cT > 1 is reached and after a 
collision the aCD molecules have no much time to es- 
cape and they rest bonded to form H-bonded aggregates 
which determine an increase of viscosity induced by the 
applied shear. This hypothesis of formation of hydro- 
gen bond clusters had been proposed in absence of shear 
as a temperature induced process [l3| . Here we propose 
that the same mechanism is activated with the applica- 
tion of shear and that the shear rate plays the same role 
of temperature in favouring the formation of H-bonded 
aggregates. 

In conclusion in this work we presented rheological 
measurements in molecular solutions undergoing inverse 
melting and we report the first experimental evidence of 
shear thickening in such a kind of systems. The viscosity 
as a function of the shear rate shows, at different concen- 
trations and temperatures, two different behaviors below 
and above a critical shear rate 7c: for 7 < 7c the fluid 
has a Newtonian behavior while for 7 > 7c it enters a 
shear thickening regime. The shear thickening transition 
is observed at critical shear rates that exhibit an Arrhe- 
nius temperature dependence jc{T) — B eyi^{Ea/ KbT) 
where Ea is the activation energy. Ea has been found 
very close to the value of the hydrogen bond energy of 
the 0-H group of the aCD molecules. This finding al- 
lows to speculate that the increase of viscosity and hence 
the shear thickening behavior derives from the forma- 
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tion of hydrogen bonded clusters favoured by the apph- 
cation of shear. This aggregation mechanism activated 
here with the appHcation of shear resembles the one sug- 
gested to explain the inverse melting phenomenon acti- 
vated by temperature (l3| . We propose therefore a single 
mechanism to be at the origin of both inverse melting in- 
duced by temperature and shear thickening induced by 
shear field in such a kind of systems. 
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